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o ■ 

' We present a general, model independent argument demonstrating that gluons produced in high 

' energy hadronic collision are necessarily correlated in rapidity and also in the emission angle. The 

, strength of the correlation depends on the process and on the structure/model of the colliding 

, ■ particles. In particular we argue that it is strongly affected (and underestimated) by factorized 

' approximations frequently used to quantify the effect. 

Q: 

ID ' I. INTRODUCTION 

^ ! 

,—1, The CMS observation of angular and long range rapidity correlations in the hadron spectrum, the so called "ridge" 
<~i ■ in proton-proton collisions [ij, has triggered a lot of discussions in recent literature 0)0- A similar if more pronounced 
Qh' correlated structure was observed in gold-gold collisions at RHIC [4] . There is a variety of candidate explanations for 
^ . the RHIC observation Q most of them utilizing strong radial flow as a collimating mechanism. Although flow 
^ ' measurements have not been reported for the LHC data, it is difficult to imagine that flow will have a significant 
. effect in p-p collisions. Thus the viable explanation should probably not appeal to any collective behavior of produced 
particles. 

The purpose of this note is to point out that at high energy, rapidity and angular correlations between produced 
^ ' particles are to be expected on very general grounds. The framework of our discussion here is similar to that of 
OO I [3|, but the argumentation will be quite general without referring to specific models of high energy evolution and/or 
0^ hadronic wave function. 

' We note that much discussion recently has been in the context of so called "glasma" 0, which refers to interacting 
, system of produced gluons immediately after the collision. We feel that placing discussion squarely in this setup is not 
' necessarily useful, and will start therefore by pointing out that the basic physics of correlations could be understood 
^""^ ] without reference to either "glasma" or "glasma flux tubes" . 

' The issue to some extent is in what frame one wishes to discuss the physics. The center of mass frame is not the 
. . [ simplest in the following sense. A very simple way of thinking about high energy scattering is the eikonal approxi- 
^ . mation, in which very energetic projectile partons scatter cikonally off the target fields. However if we are interested 
in particles produced at or close to central rapidity (like in the CMS data), we cannot use eikonal approximation in 
the center of mass frame. The incoming partons are stopped and loose most of their longitudinal momentum, thus 
' clearly the eikonal approximation breaks down. To have a valid description of particle production in the center of 
mass frame one has therefore to take into account nonlinear interactions between the produced particles, at least on 
the level of classical approximation Q. The interacting system of gluons in the center of mass frame can then be 
viewed as "glasma". On the other hand if we follow the same exact process in the lab frame, it looks quite different. 
The incoming particles are indeed very energetic and they scatter by a very small angle with p"*" ^ px- Thus recoil is 
negligible and eikonal approximation is applicable at high enough energy. Of course the physics is still the same. Thus 
for example even in the lab frame one has longitudinal electric and magnetic fields develop directly after collision, 
except that here their interpretation is rather less exotic than " glasma flux tubes" . 

The wave function of the incoming projectile carries a valence color charge density Lorentz contracted to a plane 
J~{x,x~) — p'^{x)6{x~). These charges create the Weizsacker- Williams field of softer gluons b°;{x). The WW vector 
potential is determined through the classical equations of motion[loj 

d.bUx) = p'^ix); d.b^ - d,bf - gf^'^Xix^ix) = (1.1) 

The first of these equations is essentially the Gauss' law, while the second defines & as a two dimensionally pure gauge 
vector potential. The color electric and color magnetic fields that correspond to this vector potential are localized in 
the plane containing the color charge density F^^ = bi{x)5{x~). When this configuration of charges and fields passes 
through the target, it experiences a simple eikonal rotation 

pix) ^ pix) ^ Six)pixy, b'}^b^ = Sh'^[Sp]; F+^ = b1{x)5{x-). (1.2) 
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Here iS'(a;) is a unitary matrix determined by the target color fields, which rotates the color of the incoming parton 
wave function. The effect of such a rotation is such that after the collision 



(1.3) 



Thus Gauss' law requires the presence of longitudinal electric field. In the light front the longitudinal field is not an 
independent degree of freedom, but rather is given by 
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[D,F+' - p6ix~)] = 9{x-)d, S^b,[Sp] - h[Sp 



(1.4) 



The longitudinal magnetic field is also generated by the eikonal scattering, since the vector potential after scattering 
is not a two dimensional pure gauge anymore 
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We stress again, the only physical process we are discussing is the eikonal scattering of gluons of the projectile on the 
static target. The origin of the longitudinal fields is therefore simply the Weizsacker - Williams fields of the gluons 
which are scattered out of the incoming beam. Once the direction of propagation of gluons is not parallel to the z 
axis, the color electric and magnetic fields which are transverse to the direction of propagation have a nonvanishing 
projection onto the z axis - hence the appearance of the "longitudinal fields". This is a purely kinematical effect and 
it does not involve any dynamics of formation of flux tubes, or even interaction of produced gluons as such. In the 
rest of this paper we will therefore not refer to glasma flux tubes, and our view of the scattering process will be that 
of a straightforward eikonal scattering of the projectile gluons on the color fields of the target. 



II. WHERE DO CORRELATIONS COME FROM? 



We first discuss a very simple picture of the origin of correlations, and then confront it with the available today, 
albeit incomplete formulae for the double inclusive distribution. 

Consider high energy scattering of a hadronic projectile on a stationary target in the lab frame. Since the projectile is 
very energetic, its wave function is approximately boost invariant. The boost invariance is of course only approximate, 
since at too high energy the rapidity evolution is important, and that introduces rapidity dependence inside the wave 
function. However for rapidity intervals AY < the evolution is not important^, and thus can be neglected if the 
produced particles are separated by rapidity interval which is not parametrically large. 

The boost invariance leads naturally and straightforwardly to long range rapidity correlations. Simply put, the 
incoming wave function is the same at rapidity Yi and l2- The gluon distribution at rapidity Yi and I2 are the 
same, these gluons scatter exactly on the same target, and thus whatever happens at Yi also happens at Y2. If for 
a particular target field configuration a gluon is likely to be produced at Yi at some impact parameter, a gluon is 
also likely to be produced at I2 at the same impact parameter: e voila - correlations. This is especially true in the 
context of the projectile wave function dominated by the large "classical" Weizsacker- Williams field, since in this case 
fluctuations in the wave function are small and the gluon density configuration by configuration is almost the same at 
all rapidities . This is the property of the hadronic wave function at high energy [ll[ . 

|*)=exp{z J d^xb^ix) J d7/(af (a;,r;)+<(:E,r;)) }B(a,at)|^). (2.1) 

Here ip is the wave function of valence charges, determining the distribution of the charge density p, _B is a Bogolyubov 
-type operator of the soft gluon fieds a, and the Weizsacker- Williams field b is given in terms of p by eq. (|1.1|) . For large 
projectile the WW field is parametrically large b ^ ^, while the Bogolyubov operator B produces the fluctuations of 
the gluon field of order unity. Thus for fixed p{x) the gluon density fluctuates very weakly around large average value 
determined by the classical field 

n= (fltfl) oc62 ^0( — ); (n^) - {nf ^ 1 (2.2) 

as 

The smallness of the fluctuations is clearly helpful. Although the wave function at different rapidities in a boost 
invariant projectile must be the same, the magnitude of the color field (and therefore the number of gluons) may 
differ at different values of Y for the same configuration of the valence color charge density p, if the fluctuations in 
this wave function are signiflcant. Thus in the same scattering event there may be signiflcant differences between 
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particle production at different rapidities. Still, although the quasiclassical nature of the state ea. (|2.ip ensures long 
range rapidity correlations at large values of p, it is not absolutely necessary. Even in the presence of considerable 
fluctuations in the soft gluon wave function, one nevertheless would expect positive correlations in rapidity. The only 
really necessary condition is that the density of incoming partons is large enough, so that there is large probability 
to produce more than one particle at a given impact parameter (we will quantify what we mean by "given impact 
parameter" shortly). 

Thus the long range rapidity correlations come practically for free whenever the energy is high enough so that the 
wave function of the incoming hadron is approximately boost invariant, and there is very little in the actual dynamics 
of the collision that can affect this feature. But by almost exactly the same logic we must conclude that positive 
angular correlations are also almost unavoidable. Indeed, if two gluons hit the target at the same impact parameter, 
their scattering amplitude is determined by the same configuration of the target field. Thus if the first gluon is 
likely to be scattered with momentum q, the same is true for the second gluon. One therefore expects clear forward 
correlations for gluons that scatter at the same impact parameter. Of course, the two gluons will not scatter always 
with exactly the same momentum transfer even if they hit at exactly the same impact parameter, since even a fixed 
configuration of target fields corresponds to a nontrivial probability distribution of momentum transfer. Nevertheless, 
given that this distribution has a maximum at some particular momentum transfer, the angular correlations must be 
very generic. 

The previous discussion is clearly oversimplified, since it does not address some important points. For example, 
for a soft gluon to be produced in the final state, it is not enough for it to acquire some transverse momentum. It 
also must decorrelate from the valence charge that emitted it in the incoming wave function. Otherwise it will not be 
produced as a particle in the final state, but rather as part of the Weizsacker- Williams field of the produced valence 
parton. We will therefore turn to an explicit formula that determines the gluon double inclusive spectrum in order to 
see to what extent this explicit expression is consistent with our simple discussion. 

Calculation of multigluon amplitudes at high energy has been a subject of several papers in recent years 
[i3>[i3>[i3>[3'[Ill- Nevertheless, unfortunately we still do not have a final formulae which can be used to analyse 
p-p collision at higli energy. The approach of [Ij is suited for applications when both colliding hadrons are dense, 
while that of fl^l , [l6'| , fTs*! is appropriate when one of the objects is dense and one dilute. The actual situation in p-p 
collisions at LHC is probably one where the density in the proton wave function is still not parametrically large, but 
is already not perturbatively small. In that sense the expressions of are probably more appropriate, since they do 
not discard terms which are leading at low densities, which is done in [l^. At any rate, we believe that qualitative 
features of our discussion are general enough, and we will use the explicit formulae from [iGj . 

According to (see also (isj ) the inclusive two gluon production probability is given by 

'^^ {AfHk,p)A*f{k,p))p,T (2.3) 



cPpcPkdrjd^ ' 
with the amplitude 

A^^(k,p) ^ I e^'^^+^f"^ [u{z-x,)[S{x,)-S{z)rp^{x,)][f,{u-x^)[S{u)-S{x2t%\x2)] 
Mz - xi)fj{u - xi){ [S{xi) - S{z)] p{xi) [SHu) + S^ixi)] y 

1 

Mz-u)f,{u-x,){{S{z)-S{u))p{x,)SHu)}'"' . (2.4) 
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Here 



and we have defined p = —iT°'p°^. The charge density is normalized such that for a single gluon p° = gT"". In these 
formulae p°'{x) is the valence color charge density in the projectile wave function, while S°'^{x) is the eikonal scattering 
matrix determined by the target color fields. The average in ea. (|2.3p denotes averaging over the projectile and the 
target wave functions. We also note that in this expression the gluon with momentum p is assumed to have larger 
rapidity, and thus the emission of the two gluons is not completely symmetric. 

The actual eq uation given in (l6| is slightly different and this requires explanation. The color charge density 
appearing in |16l| are quantum operators with quantum commutation relations of SU{N) algebra. On the other hand 

in ea. (|2.4p are c-number functions, and averaging over p is understood as averaging over a classical ensemble 

{0)p = [ DpWp[p]0 . (2.6) 
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As discussed at length in [T^l, the full quantum averaging is equivalent to classical averaging procedure ea. (l2.6p for 
totally symmetrized products of operators /5° . To obtain eq. (|2.4p from the expression given in we have rewritten 
the original expression of [l6| in terms of the anticommutator of two p's in the first line. Thus the "classical" expression 
p^p^ is equivalent to the quantum operator ^{/5'*,/5^}. This procedure of reordering generates an additional term in 
ea. (l2.4p which is not present in [l6| . 

The physical meaning of the three terms in ea. (j2.4p is straightforward. The first term corresponds to independent 
production of the two gluons. This term is leading in the limit of large color density p ~ I/17. One should keep in 
mind, however that in this limit other terms not included in eq. (|2.4p are equally important The second term 

corresponds to production of two gluons emitted from the same color source in the incoming projectile wave-function. 
The third term corresponds to the process whereby the softer gluon has been emitted in the wave function by the 
harder one, with both gluons subsequently produced in the collision. In terms of BFKL ladders, the (square of the) 
first term is a part of the diagram containing two independent ladders, while the (square of the) last two terms 
describe emission of two gluons contained in the same BFKL ladder. 

To calculate the cross section one has to square the amplitude. This produces many terms, since in accordance with 
our previous discussion one has to symmetrize the factors of charge density between the amplitude and the conjugate 
amplitude. The full expression for the cross section is given in the Appendix A. Here in the text we only reproduce 
one part of this expression which arises from squaring the first term in the amplitude eq. (j2.4p which is responsible for 
independent production of the two gluons. 

'^^ (a4+a2+a3)p,T (2.7) 



with 

= / e^M-«)+*P(«-^) /■ /(z-xi)-/(xi-z)/(M-i2)-/>2-u) (2.8) 

X {p{xi)[S\xi) ~ S\z)]\S(xi) - S{z)]p{x^)] {p{x2)[SHu) - SHx2)][S{u) - S{x2)p{x2)} 

The explicit expressions for and are given in the Appendix A. Although these are long expressions, only 
the term given in eq. (j2.8p is relevant for most of our discussion. This is precisely the contribution which has the 
physics discussed above. The two gluons here are produced independently from each other, but from exactly the same 
configuration of sources through scattering on the same target field. In the other terms the two gluons are correlated 
with each other in the incoming wave function, and thus these terms contain additional physics. Let us therefore 
consider the cross section of eq. (|2.8p . It is very easy to see that it indeed produces angular correlations. One can 
write it as 

<j^{k,p) = {<j{k)a{p))p,T (2.9) 

where 

<T{k) = / e"=(^-^") / f{z - xi) ■ fix^ - z) {p{x^)[S\x,) - S\z)][S{xi) - S{-z)]p(x^)} . (2.10) 

For fixed configuration of the projectile sources p{x) and target fields S{x) , the function a{k) as a function of momentum 
has a maximum at some value k = q. Therefore clearly the product in eq. (|2.9p is maximal for k = p = q. The value 
of the vector q of course differs from one configuration to another, but the fact that momenta k and p are parallel 
does not. Therefore after averaging over the ensemble CT^(fc,p) has maximum at relative zero angle between the two 
momenta. 

We reiterate, that even though averaged over all configuration {a{k))p,T must be isotropic, there is absolutely no 
reason for it to be isotropic for any given configuration. The strength of the maximum of course depends on the 
detailed nature of the field configurations constituting the two ensembles (the projectile and the target). We will 
discuss some qualitative features of these in the next section. But first, it is interesting to ask is the maximum of cr(fc) 
unique, or perhaps there is finite degeneracy. It is in fact easy to see that the maximum is doubly degenerate. The 
probability cr(fc) can be written in terms of the single gluon production amplitude a{k) 

a1{k)= [ dze^'^^ [ f,{z-x{)[S{xi)-S{z)f''p\xi). (2.11) 

J J x\ 

Since the amplitude a is real in coordinate space, we have 



(7{k) = a{k)a*{k) = a{k)a{~k) . 



(2.12) 
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Configuration by configuration tfiis is clearly symmetric 

cr(fc)=CT(-fc). (2.13) 

The " classical" ' contribution to the two particle inclusive production probability is therefore symmetric under 

o\k,p) = a\-k,p) . (2.14) 

and must have two degenerate maxima - at relative angles A(/) = 0,7r. The feature of the equal strength of forward 
and backward correlations stays true beyond the specific form of the single gluon production amplitude ea. (|2.11l) . 
In the limit of dense projectile eq. (|2.10|) is not strictly applicable. However, although one does not have an analytic 
formula valid in this limit, the formalism of [T4| expresses the double gluon production probability in terms of a 
solution of classical equations of motion. This solution, which is precisely the single gluon production amplitude, is 
real in coordinate space for any fixed configuration of projectile and target fields. The reality of the amplitude in 
the coordinate space is the only property required to establish degeneracy, and thus degeneracy persists in the dense 
projectile limit as well. We note that numerical results of [3] do not seem to display the exact symmetry of cq. p.l4p . 



which may be an indication of some subtlety of the numerical procedure of [18 1. 

The other terms in the gluon production amplitude ea. (j2.4[) are not likely to lead to correlations of the type just 
described. 

The third term in eq. (j2.4p . where the gluon produced at the point z is emitted from the other observed gluon at 
the point u, disfavors production at the same impact parameter because of the suppression factor S{u) — S{z). The 
two gluons when scattered at the same impact parameter do not decohere, but rather scatter as a single coherent 
state, with the gluon at z emerging in the final state as part of the Weizsacker- Williams field of the gluon at u. On 
the other hand whenever the two gluons do decohere, since they were correlated in the incoming wave function, they 
emerge in the final state with large relative transverse momentum. Thus this particular term in the amplitude mostly 
leads to back to back production in the final state and is responsible for the large away side, rapidity independent 
maximum at relative angle tt, prominently present in the data. 

The second term in ea. (|2.4p favors production at the point u close to xi, but z far from xi. Thus one expects the 
momentum of the gluon produced at z to be uncorrelated with that of the gluon produced at u. Whenever the gluon 
at u is produced with significant transverse momentum, the balancing transverse momentum resides at the "valence" 
rapidity. This term is therefore responsible for the away side peak between one of the observed particles and another 
particle produced at a more forward rapidity. 

One can estimate the overall magnitude of the correlation by the following simple argument. In order for two 
produced gluons to be correlated in the final state, they have to be close in the initial state and also scatter off the 
same target field. We will assume that both the target and the projectile are characterized by corresponding saturation 
momenta Qs ■ The inverse of the correlation momentum is the correlation length in the hadron L ^ 1/Qs- It is 
reasonable to expect that typical field configurations contributing to the hadronic ensemble of, say the target, have 
variation only on distance scale greater than 1/Qs- Thus the two gluons that hit the target at distance x < 1/Q^ 
apart from each other scatter on the same field. By the same argument, for the two incoming gluons to be in the 
same state they have to be located in the impact parameter plane no further than 1/Qf away from each other. Thus 
for correlated production the two gluons need to be within the radius q^^^ of each other, where Q™"^ is the larger 

of the two saturation momenta Qf and Q^. On the other hand the total number of produced gluons is proportional 
to the total transverse area of the smaller between the two objects participating in collision. Thus parametrically 



d^N dN dN 



d^pd^k d^k d^p 



dN dN 1 



This estimate is of course parametrically the same as given in 

We would like at this point to make contact with the recent paper ^ . The calculation of gluon production in fsj 
is based on simplified version of eq. (|2.9p supplemented with specific prescription for averaging over the projectile and 
target fields. Specifically, [3[ expands the scattering matrix S to first order in target fields, and keeps only the leading 
term S{x) — > 1 + a{x). The expression for then becomes a homogeneous function of the target and projectile fields 

cr'* ^ [paap){paap) (2-16) 

For simplicity we suppress the color indices and transverse coordinates on all the functions. One next averages over 
the charge densities assuming Gaussian ensemble 

(P^)=3(p2)(p2) (2.17) 
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and similarly for the target. And finally the high energy evolution is included by substitution 

{p{x)p{y))^^{x-y) (2.18) 

with $ taken to be a solution of the Balitsky-Kovchegov equation [l9|. Although the angular distribution has not 
actually been calculated in the authors argued that the correlation should in fact have a maximum at collinear 
momenta. The subsequent numerical investigations based on the same approach bear this out [20j . 

Our general discussion provides an intuitive explanation for this result and also makes it clear that the presence of 
the correlations does not depend on the specifics of the approximation used to estimate Jif^^-ip ■ The magnitude of the 
effect however may depend on the approximation quite strongly. We next want to comment on this issue. 



III. ISSUES WITH AVERAGING. 



From eas. (|2.9l2.10p we know that the basic averages that one needs to calculate are of the type 

{[S\x)S{z)r'[S\y)S{u)r% (3.1) 

and similarly for the projectile 

{p'^{x)p\x)p^{y)p\y))p. (3.2) 

The Gaussian averaging procedure described above is fairly restrictive, in the sense that as any Gaussian averaging it 
probably tends to underestimate correlations. In particular Gaussian averaging over color singlet ensemble necessarily 
puts the densities in eq. (|3.2[) pairwise into color singlet states. As pointed out in [2l[ this leaves out some possible 
configurations which are overall color singlets, but where no two factors of p form a color singlet separately. This for 
example happens, when the factors of p are pairwise in color octets, with the two octets forming an overall singlet. 
Such configurations in principle can also contribute to the correlated part of the particle production. Formally they 
are suppressed in the large N,. limit. However the correlated part of the production probability itself when calculated 
with the Gaussian averaging is also suppressed by l/A^^ relative to the uncorrelated part, and thus omission of these 
terms may be dangerous [21'|. Physically these terms correspond to interference contributions. For example when the 
two factors of p in ea. (l2.10p are in an octet, this corresponds to a situation when the charge densities in the amplitude 
and complex conjugate amplitude are different, but still the same gluon in the final state is produced due to the 
difference in the scattering factors S in the amplitude and the conjugate amplitude. 

Although these l/jV^ suppressed terms are interesting, taking them properly into account requires one to go beyond 
the dipole model [2^ and the BK equation, and in the dense region studying the full B-JIMWLK evolution (23| . 
However it is not obvious that even in the leading order in 1/Nc the Gaussian approximation is adequate to discuss 
correlated production. Here we would like to discuss only these leading order terms. We will argue that Gaussian 
averaging procedure is likely to miss terms in the correlated production probability which are of the same order in 
1 /Nc as the uncorrelated piece. 

The leading Nc piece in ea. (l2.9p comes from the configuration where the charge densities in each one of the single 
gluon production amplitudes are in the color singlet. The relevant average to calculate is 

{p'^{x^)p^{x^)p'{x^)p\x2))p{^r{[SHx^) - S\z)][S{x^) ~ Siz)]) T,{\S\x^) S\u)\\S{x-,) ~ S{u)W)t . (3.3) 

Let us first concentrate on the projectile average. As mentioned above, averaging with a Gaussian weight one obtains 
in the leading order in 1 jN^ 

{p''{xi)p''{xi)p''{x2)p''ix2)) Gaussian = {p" (xi) (xi)) Gausstanip'' {X2) p'' {X2)) Gaussian ■ (3.4) 

In this approximation therefore, clearly the correlated piece in the production probability vanishes, and only the 
subleading in 1/Nc correction resurrects the correlations. We stress however, that this is not the result of the leading 
Nc approximation per se, but rather of the Gaussian averaging procedure. 

It may be tempting to think that factorization in the large Nc limit is natural due to presence of large number of 
degrees of freedom, and therefore in some sense large Nc might act similarly to heavy nucleus. However this is not 
the case. Even though the number of degrees of freedom is large, even in the large Nc limit the theory has legitimate 
states which contain small number of particles. A color dipole is an example of such state. It is a superposition of 
many states (different color orientations) of two particles, rather than a state with many particles. In a state like this 
the central limit theorem does not hold, the fluctuations in density can be large even in the large Nc limit, and it is 
the large fluctuations in the ensemble that break factorization of correlation functions. 
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In fact a very similar question was considered a while ago in [24j in connection with factorization of dipole densities 
in the dipole model 22]. Indeed the observable we are interested in eq. (|3.3p is rather similar to the dipole density 

n(a:i,xi) - (3.5) 

As shown in [24| within the dipole model (which is defined entirely within the large Nc limit [1^) the product of two 
densities does not factorize, but rather behaves as 

{n{xi,Xi)n{x2,X2)) - {n{xi,X2)){n{x2,X2)) ^ {n{xi,X2)){n[x2,X2))b^^ (3.6) 

where h is the transverse distance between the two dipoles and A is a number, whose exact value is unimportant for 
us. This result hods in the limit where the distance between the two dipoles is much greater than their respective 
sizes, and thus it does not display any angular correlation between the orientations of the two dipoles. Nevertheless 
ea. (l3.6p clearly exhibits the fact that factorization is not an inherent property of the large Nf. limit. Once we accept 
that the factorization is broken, it is natural to expect that the actual correlation function in the regime where the 
two dipoles overlap in space, also exhibits angular correlations in the orientation of the two dipoles. 

Note that this is precisely the regime relevant to our discussion of angular correlations in emission. The same 
configuration of color charges produces the same gluons (at different rapidities) , which produce correlated hadrons in 
the final state. Thus the most important region of the phase space is when all four points in the correlator eq. p.3p 
are close to each other, in the sense that they are all within the correlation length 1/Qs- It is very hard to imagine 
that in this regime factorization holds (see [2^ [2^ for more discussion of such correlations) . Thus we indeed expect 
that any realistic non-Gaussian weight function for the ensemble averaging will lead to a nonvanishing contribution 
to the correlated piece of gluon production even in the large Nc limit. 

Turning to the target averaging in eQ. (|3.3l) . the terms that have to be averaged are of the type of observables 
described in the large Nc limit by the dipole model [1^ 

(Tr { [^t ^^)s{z)] } Tr { [S^ {y)S{u)]}}T = {s{x, z)s{z, x)siy, u)s{u, y))T (3.7) 

where s{x,y) — Tr[Sp{x)SF{y)] - is the scattering amplitude of the fundamental dipole, and the equality in ea. p.7p 
holds in the large Nc limit. The approximation which is frequently used in the literature to calculate the averages of 
this type also invokes factorization 

{s{x, y)s{u, v)) = {s{x, y)){s{u, v)) . (3.8) 

The target averaging of W\ would follow from this approximation in the limit of weak fields. When the dipole {x, y) 
is far from the dipole {u,v) (much further than 1/Qs), the factorization is a good approximation since the fields on 
which the scattering amplitude is calculated are not correlated with each other. However, as before we are clearly 
interested in the case where all the points are within the distance of order 1/Qs or smaller. In this case, just like for 
the projectile, the factorization of ea. (|3.8p is not a property of the large Nc limit but is rather an ad hoc assumption, 
used only due to its simplicity. 

Strict factorization of the type eq. (|3.8p is only possible if the statistical ensemble consists of a single configuration. 
There is however no reason to expect that in the large Nc limit fluctuations around some leading configurations are 
suppressed by powers of l/Nc- For example, the energy evolution of s is given by the dipole model Hamiltonian, which 
does not contain Nc at all. The probability distribution of the dipole model W[s] evolves with rapidity according to 

[Ulli] 

with a - the 'tHooft coupling, which is finite at infinite Nc- Thus any nontrivial initial distribution W[s] evolves 
smoothly to higher energy and remains nontrivial. 

In fact it would be very interesting to see what happens to correlations that one can put into the initial ensemble, as 
rapidity grows. Technically one should choose an ensemble Wq[s] of initial configurations s{x, y), which contains short 
range correlations. These correlations should be confined to within the saturation radius, that is {s{x,y)s{u,v)) ^ 
{s{x, y)){s{u, v)) only when |a; — u|, |?/ — w| < 1/Qs- Each configuration of the ensemble should be evolved independently 
according to the BK equation [2^ . The correlations at the final rapidity are then calculated by averaging the correlator 
calculated in the final ensemble over the ensemble of initial conditions. 

This procedure is the same as the one implemented in (26| . However in order to study the angular correlations 
between the produced particles, it is important to understand angular correlations in the target ensemble. Therefore 
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the individual members s{x — y) of the initial ensemble should not be isotropic. The rotational invariance must be 
restored by averaging over the whole ensemble rather than configuration by configuration. 

As the rapidity grows, so does Qs and the correlation radius shrinks so that one naively expects only the points 
inside this new saturation radius l/Qs{Y) to stay correlated. We note however that it is also possible that the 
correlation length in this particular channel is larger than 1/Qs- For example it was shown recently in [27] in a QCD 
like model, that the saturation momentum itself varies rather slowly in the impact parameter plane. In some sense one 
can ascribe to it a correlation length which is far greater than 1/Qs, and is instead A oc ^e'^'" Although there 

is no particular reason to expect that such a large correlation length will also determine the scale of disappearance of 
correlations of the dipole amplitudes, if it does it would be very interesting and one would have to revise the estimate 
ea. (|2.15p accordingly. 

To summarize, there are good reasons to expect that the factorization of both projectile and target averages is 
broken at leading order in 1/Nc in the kinematical domain relevant to the correlated production of particles. To 
study this question one certainly has to go beyond simple rotationally invariant solutions of the BK equation. While 
technically challenging, it would be very interesting to understand and quantify this effect. 

IV. CONCLUSIONS 

We conclude with a short summary. We have shown that rapidity and angular correlations are a very general feature 
of particle production at high energy. They are an automatic consequence of the boost invariance of the projectile 
wave function, provided two conditions are met. 

First, there should be some "classical" component in the (projectile) wave function, meaning that the short range 
rapidity fluctuations in the wave function should not be overwhelming and should not lead to complete decorrelation 
between gluons at different rapidities. This condition is satisfied in QCD for dense as well as dilute projectiles. In 
the dense case the soft gluon wave function is dominated by the " classical" rapidity independent Weizsacker- Williams 
field. For dilute projectile fluctuations are significant, but their contribution to gluon density is not parametrically 
larger than that of the WW field. 

Second, the particle density/scattering probability must be large enough so that there is non-negligible probability 
to produce more than one particle at fixed impact parameter within the saturation radius 1/Qs- This condition is 
satisfied when both colliding objects are dense and the values of the two saturation momenta are not vastly different. 
In that case the number of gluons in the incoming wave function within area l/Q^ is of order l/as{Qs), while the 
probability of producing one particle is of order one. For scattering of two dilute objects this condition is not met. 
However at high enough energies the density in the wave function grows and the effect must become visible. 

The angular correlations originate from configuration by configuration fluctuations of the projectile and target 
structure in the transverse plane away from a rotationally invariant state. The effect we have discussed here has 
several tell-tale features. Produced particles are correlated in angle, with forward and backward correlations being 
of equal strength in the case where the two colliding objects are nuclei. When the colliding objects are not dense, 
there is an additional contribution to particle production, from a " single ladder" which significantly enhances back to 
back correlations. This contribution is responsible for the bulk of the observed back to back correlated production. 
The correlation is present also in the magnitude of the transverse momentum and not just in the angle. The latter 
correlation in fact does not require local rotational asymmetry of the projectile/target configurations. It would be 
interesting to try and measure these correlations as well. 

The relative magnitude of the forward correlations should initially increase with energy for p-p collisions, since the 
relative importance of the single ladder terms diminishes. Interestingly however, the estimate of eq. (j2.15p suggests 
that the effect decreases with energy, once the colliding systems can be treated as saturated objects with well defined 
saturation momentum, since the saturation momentum grows with energy. Thus at very high energies the effect should 
disappear. If we apply this logic also to nuclear collisions, we should conclude that the effect if observed by ALICE 
should be significantly smaller than that observed by PHOBOS and STAR. Our discussion of course disregards the 
effects of flow, which are generally believed to be very important for nucleus-nucleus collisions. The latest STAR data 
[soj support this view. It is possible therefore that our considerations about angular correlations are not valid for 
nuclear collisions, in the sense that the main mechanism of collimation is indeed due to the flow. It would nevertheless 
be interesting to try and disentangle the flow effects from the intrinsic correlations in the initial state discussed in 
this note. We also note that the estimate ea. (|2.15l) refers not only to angular correlation, but rapidity correlation 
in general. Thus independently of the question of radial flow, if the observed long range rapidity correlations are 
due production from correlated domains in the boost invariant incoming wave function, the trend should be that of 
decreasing correlated production going to higher energy. 

Finally, we have argued that the correlations must survive also in the leading order in expansion. Their 

subleading nature in current numerical implementations is due to the factorization assumption which is not valid in 
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the region of the phase space relevant for the correlated production. We believe that improvement of this aspect of 
current calculations is imperative in order for the results to be quantitatively reliable. 



Appendix A: Ordering 



In this Appendix we derive the completely symmetrized expression for the doubly inclusive gluon production. When 
averaging this expression over the projectile wave function, the color charge density p should be treated as classical 
commuting variable. Simply squaring the expression for the amplitude in ea. (|2.4l) . one gets eq (|2.7p with cr'* given in 
eq.dlH) and 



^ik(z — z) -\-ip{u—u) 



u,z,u,z 
2 



(Al) 



— /(z - Xl) ■ f{xi - z) f{u ~ Xl) ■ f{xi - u) 



xTr{ [S{u) + S{xi)] p{xi) [S\xi) - S^z)] [S{x,) - S{z)] p{xi) [S^{u) + S\xi)] } 



+9^f{z - u) ■ f{u - z) f{u - Xl) ■ fixi - u) 



xTr 
2 

xTr 



{S{u)pixi) (SHz) - S\u))} {{S{-z) S{u)) pi^i)SHu)} 

f{z - Xl) ■ f{u - z) f{u - Xl) ■ f{xi - u) 

{S{u)p{xi) (SHz) - S\u))} [S{xi) ~ S{z)] p{xi) [S^{u) + S^xi)] 

f{z - U) ■ f{xi - Z) f{u - Xl) ■ f{xi - u) 



xTr 



[Siu) + Sixi)] pixi) [S\xi) + SHz)] {iS{z) - Siu)) pixi)SHu)} 



^3 _ ^ik{z-z)+ip{u-u) / ^^2) 

-|/(z - Xl) ■ f{xi - z) f{u - Xl) ■ f{x2 ~ U) 

xp{xi)[SHxi) ~ S\z)\ \S{xi) - S[z)\ pixi) {S^{u) + S\xi)\ \S(u) - S{x-,)\p{x.,) 
-|/(^ - Si) • - z) f{u - X2) ■ f{xi- U) 

xp{x2)[S\u) - SHx2)] [S{u) + S{xi)] p{xi) [S\xi) ~ S^{z)] [S{xi) - S{z)]p{xi) 
+gf{z - u) ■ f{xi - z) f(u -Xl) ■ f{x2 - u) 

xpixi)[S^ixi) - S\z)\ {(5(z) - S{u))p{xi)S\u)) {S(u) - S[x2)\p{x2) 
+gf{z ~ Xl) ■ f{xi - z) f{u - X2) ■ f{xi ~ u) 

xpix2)[S^iu) ~ S\x2)\S{u)p{xi) {S\z) - S\u)\ {S{xi) - S(z)\p(xi) . 

However in this expression the factors of p coming from the amplitude are not symmetrized with respect to factors 
coming from the conjugate amplitude. To derive the fully symmetric expression we instead follow the formalism of 

The amplitude of a single gluon production in the dense-dilute scattering is given by 

Q1{z) = J d'xMz - x) [S'^'iz) ~ S''\x)] 4[S,x] , r^{S,x) = -tr !^S{x)T'^j^^ (A3) 

In this expression the projectile wave function has been "integrated out" [13]. It can be however "integrated back in" 
by representing the right rotation operatorsd as the projectile charge density operators acting on the projectile wave 
function 25[ : 

JbMP) = Px\P), rn{x)Jl,{y)\P) = pIpI\P) (A4) 
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The operator associated with inclusive production of a single gluon with momentum k is 



i k (z — z) r^a 



2ti 2ti 



(A5) 



where S is the single gluon S-matrix in the amplitude and S in the conjugate amplitude. Double gluon production 
without rapidity evolution between the two produced gluons is given by 



dN 



cPpcPkdrjd^ 



{Og{k)Og{p))p.T\s = S 



(A6) 



The operators Qf{z, [S]) do not commute with each other. We will move all the right rotation operators Jn to the 
right and then convert them into the operators p according to eq. (jA4[) . We also keep all operators Jr[S] to the right 

0fJR[S]. 



The amplitude A^j" reads 



^ikz+ipu f^a / 



{z)Q){u)^ / e»'=-+'f" / [h{z-x^)[S,-S^,Y-f,{u-X2)[Su-S.,,]'"'j'k{xi)Ji{x2) 



J Z^U 'J Z^U J XlX2 

+ h{z - x^)[S. - S,,r - X2)[Su6u., - T;^,p^^ (A7) 

Symmetrizing the operators within this amplitude, one recovers eq. (|2.4p . Our aim here is to obtain the fully ordered 
expression for the probability. We therefore square the amplitude first and then perform the full symmetrization 
between the operators p entering both the amplitude and its conjugate. Let us introduce some notations that will 
help to make our expressions compact. 



f{z - Xi) ■ f{xi - z) 



(A8) 



^zxizxi — ^zxizxi [{Sl Sl^){Sz Sxi)] 



^XlXl (^) 



^ik{z — z) zpab 



= Fl'^.A-k) (A9) 



Squaring the amplitude eq. (jA7l) we obtain 

dN 



ik{z-z) riab njab — M rotCC_„C'_ \]ab 

^zxizxi J ^"zxizxi — ^^zxizxi [^z \'-'z '-'xiJl 

(AlO) 



d^pd^kdrjd^ 

with (sum over both the color and coordinate indices is implied) 

= Ff,,Ak)Fxt,Ap)PxJxjijk 

y^3 _ / ip(u-u) pab (U\ ]\/T^d __ rpa _ pab ( U] nf^d (^\rpa 

7 J- uXl \ / UX2UX2 cm XlXl \ ) ^X2X2 V-r^/ CT 



(All) 



'■c p^d pb 

rX2 rX2 rXl 



U 

^ip{u F^^.-^^{k) Nux2UX2^c'm [i^u^uxi ~ Sl^Sx-i^x2){SuSuxi —3x2^X1X2)] ^nc/ ^^^2 (-^12) 



Following |17'], the operators are substituted by commuting c-number functions according to: 

-1 ba 



Pxi Pxi 



5p%, J 12 V 5pi, 



(A13) 
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After long algebra we obtain the fully ordered expression 



'^^ ={a^ + a^ + a')p,T (A14) 



cPpcPkdrjd^ 
with 

a" = [pFik)p] [pF{p)p] (A15) 

= + g{^PxFxy{~k)pyFyy,{-p)py, + ^PxFxy{k)pyFy^{p)p.u] + ^tr[pxFxx{k)][pyFy^{p)p.u]] + 
+ ^[PxFxy{k)py] tr[pxFxx{p)]^ 

= 9(^-^PxFxy{k)pyFyw{-p)puj + -PxFxy{~k)PyFyui{p)Pw + P X Fxy {- k) PyG ym {p) Pw " P X Fxy {k) fiy G ym {- p) Pw " 

g«p(«-a) [p^^Fx,u{~k)px^Mux2ux2Px2 + Px ^ Fx ,u (k) Px ^ M^x 2UX2 PX2]) (A16) 



= a^+g^[- -tr[pxFxxik)T''][Fxy{~p)pyr ~ ^tr[FxxXk)T''][pxFxy{-p)pyr - j^tr[T''Fxx{k)px][Fxy{-p)pyr 

+ ^[PxFxy{k)Y'tr[PyFyyip)T'] - ^[p x Fx y {- k)]'' tV [fjy Fyy {p)T + ^tr[p x Fxy (k) f) y Fy x (p)] + 

+ ^tr[Fxx{k)T''][pxFxy{p)py]'' + ^tr[pxFxx{k)]tr[pyFyyip)] + ^[pxFxy{-k)py]''tr[T'Fyy{p)] - 

- ^[PxFxy{~k)]-tr[T-Fyy{p)py] + ^ i r [^^ F,, ( fc ) [F,^ (p) ] + ^tr[T«F,,(fc)p,] - 

- ^[PxFxy{k)rtr[T-Fyyip)py] + ^[p X Fxy (k) fly^ triT^ Fyy (p)] ~ ^ip X Fxy k)]- tV [T'^ G yy (p) fJy] + 

+ hr[T''Fxxik)px][Gxyip)pyr - l[PxFxy{k)]''tr[T''Gyy{-p)py] + ^tr[T''Fxxik)px][Gxyi~p)pyr + 
\ fe'P^^-^^ {[p,,F,,u{~k)rtr[T-Mux2Ux2Px2] - ir-[T'^F„,,(fc)p,J[M„,,a,,psJ'^ - 



UU 



2 

tr[px,Fx,nik)T''][Mux,ux2Px2V + [Px,Fx,u{kTtr[T'"Mux2ux2Px2 



2 2 



^tr[pxFxy{k)pyFyx{-p)\ - ^tr[Fxy{-k)pyGyx{p)px] + ^tr[Fxy{k)pyGyx{-p)px 



I G ^ (tT\^Px2 Fx^U ( ^^Px2 ■^UX2UX2 ] ^^[Px2 -^3^2 u{}^')Px2 ■^'^UX2UX2 ] ) 



2 



^ c\ I ^ { -^UX2UX2^'^[Fuu{ k^ Px2^u^uPx2] ^UX 2UX2^^[Fx 2u{_ Px2 X2 PX2] 

^ J UU ^ 

~ NuX2UX2tf\FuX2i. k)Px2S'x2^'0-Px2\ ^" A^^jj^jfl [-Pkl 1:2 ( ^ ^)P2;2 '^£!;2 "^^S ] ^ (-'^l''') 

This can now be averaged over the projectile wave function with the help of classical ensemble averaging procedure 

(...)P = / [dp]Wp[p]... (A18) 
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